Dense Crystalline Dimer Packings of Regular Tetrahedra 



Elizabeth R. Cheif] 



Department of Mathematics, University of Michigan, Ann Arbor Michigan 48109, USA 

Michael Enge|^ and Sharon C. Glotzeij^ 

Department of Chemical Engineering, University of Michigan, Ann Arbor Michigan 48109, USA 
(Dated: Submitted: December 27, 2009; Revised: July 27, 2010) 

We present tlie densest known packing of regular tetrahedra with density (j) = 
= 0.856347 . . . Like the recently discovered packings of Kallus et al. and 
Torquato-Jiao, our packing is crystalline with a unit cell of four tetrahedra forming 
two triangular dipyramids (dimer clusters). We show that our packing has maximal 
density within a three-parameter family of dimer packings. Numerical compressions 
starting from random configurations suggest that the packing may be optimal at 
least for small cells with up to 16 tetrahedra and periodic boundaries. 
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I. INTRODUCTION 



The problem of the packing of tetrahedra, which has modern-day applications to such 
wide-ranging topics as metamaterials with novel optical properties, nanomaterials [IH3], and 
virus formation, hails back to the early Greeks [1]. Aristotle, in discussing the assignment 
of geometrical figures to "heavenly bodies", mistakenly believed that regular tetrahedra 
("pyramids") tile (Euclidean) space perfectly. Nearly 1800 years later, Johannes Miiller 
(aka Regiomontanus, 1436-1476) contradicted Aristotle's claim [5l|6]. 400 years after Miiller, 
Minkowski showed that the densest lattice packing of any convex body must satisfy certain 
constraints [7] and mistakenly argued that the densest lattice packing of tetrahedra (i.e. one 
tetrahedron per unit cell) had density = ^ [H]- In 1900, Hilbert posed the problem of the 
densest lattice packing of regular tetrahedra as a special case of the 18th of his famous list 
of problems |9]. Gromer conjectured ^0\, and Hoylman later proved [H], that the densest 
lattice packing of a single tetrahedron is = ^. 

In 1972, Stanislaw Ulam conjectured that spheres would have the lowest maximum pack- 
ing density of all convex bodies, including tetrahedra [12]. That density is ^sphere = = 
0.740480... [I3|. In 2000, Betke and Henk developed an efficient computer algorithm to 
compute the densest lattice packing of any convex body, and applied it to the Archimedean 
solids [T^. Conway and Torquato [T5] used Betke and Henk's algorithm to compute the 
packing density of tetrahedra derived from the densest lattice packing of icosahedra, and 
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also examined other promising packings. All of them packed worse than spheres, with a 
maximum packing fraction ^ 0.7175. This led them to suggest the tantalizing possibility 
that the Ulam conjecture might, in fact, fail for tetrahedra. They further proposed that 
tetrahedra might have the lowest packing density of all convex bodies. Motivated by this, 
Chaikin and coworkers [12] performed experiments on tetrahedral dice packed in spheri- 
cal and cylindrical containers. They reported random packings with densities as high as 
0.75 ± 0.03 [16] and 0.76 ± 0.02 [17], with the error arising from the rounded dice vertices 
and edges [T7] . 

In 2008, Chen proposed a packing of nonamers (nine tetrahedra forming two pentago- 
nal dipyramids that share one tetrahedron) arranged in layers, whose density (0 ~ 0.7786) 
clearly exceeded, for the first time, the maximum packing density of spheres p!8j. Since 
Chen's publication, a flurry of activity has resulted in successively higher and higher pack- 
ing densities for tetrahedra. It was soon demonstrated that the nonamer crystal could be 
numerically compressed to slightly higher packing densities (0 = 0.7837 [I9] and = 0.7820 
[20]). The first disordered (i.e. non-periodic) packing of tetrahedra to exceed the maximum 
packing density of spheres was reported by Haji-Akbari et al. [T9l 121] : they obtained a 
packing density of = 0.7858 using Monte Carlo simulation of systems of 8000 tetrahedra 
compressed from random initial conditions. Their disordered packing contains a preponder- 
ance of pentagonal dipyramids, arranged randomly and thus differently than in the Chen 
structure, along with other motifs such as icosahedra. 

Previously, all dense ordered packings proposed for tetrahedra were based on analytical 
construction (or numerical compressions thereof). Haji-Akbari et al. showed, using MC 
simulations of initially random systems containing up to nearly 22,000 tetrahedra, that at 
packing densities > 0.5 an equilibrium fluid of hard tetrahedra spontaneously transforms 
to a dodecagonal quasicrystal, which can be compressed to = 0.8324 [21]. By numerically 
constructing and then compressing four unit cells of a periodic quasicrystal approximant 
with an 82-tetrahedron unit cell, they obtained a packing density as high as = 0.8503. 

Motivated by a numerical search, Kallus, Elser and Gravel [22] found a one-parameter 
family of dimer cluster packings obtained via an analytical construction with a density of 
= = 0.854700 . . ., exceeding the density of the quasicrystal approximant. Each dimer 
cluster contains two face-sharing tetrahedra, and is equivalent to a triangular dipyramid 
(or bipyramid). Two dimers comprise a single unit cell. Torquato and Jiao generalized the 
analytical construction of Kallus et al. to a iwo-parameter family of packings, and obtained 
an even denser packing with = = 0.855506 . . . |23j. 

In this paper, we further generalize the Kallus et al. and Torquato- Jiao family of dimer 
cluster packings to consider a i/iree-parameter family of packings. We obtain a maximum 
packing density of = |pj = 0.856347 ... as the optimal solution within the three-parameter 
family and show that the Kallus et al. and Torquato- Jiao packings are special cases of our 
more general construction. Furthermore, we present numerical compression simulations of 
small systems of tetrahedra using standard isobaric Monte Carlo with variable simulation 
cell shape. For systems with 4, 8, 12 and 16 tetrahedra with initially random configurations, 
the analytically predicted dimer crystal with = 0.856347 is recovered. 
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II. ANALYTICAL CONSTRUCTION 

A. Double dimer configurations 

Definition 1. A dimer (cluster) is a basic building block formed by two face-sharing regular 
tetrahedra arranged into a triangular dipyramid. The vertices of a positive dimer are chosen 
up to translation as 

o=(+2,+2,+2), 

p = (+2,-l,-l), g = (-l,+2,-l), r = (-l,-l,+2), (1) 
s = (-2,-2,-2), 

such that p, q, r span the common face. A negative dimer is related to a positive dimer by 
inversion. The positive dimer with vertices a, p, q, r, s centered at the origin is called +F2; 
the corresponding negative dimer with vertices —a, —p, —q, —r, —s centered at the origin is 
called — F2. 

Definition 2. A crystalline arrangement of tetrahedra with four tetrahedra in the unit cell 
forming one negative dimer and one positive dimer is called a double dimer configuration. 
// the configuration is free of overlaps, then it is a double dimer packing. 

The set of possible double dimer configurations P is a set of translates of the two dimers. 
We place the positive dimer +F2 at the origin and translate it by the lattice L+. The 
negative dimer — F2 is placed by the lattice coset L^. Without loss of gencrahty we can 
assume that L+ is spanned hy a + b, b + c, c + a and choose the offset of L~ as d + a 
(equivalently d + b, or d + c), where the lattice vectors (a, b, c) and the offset vector d are 
free to vary. This means that 

L"*" = {riaCL + Ubb + UcC \na + nb + nc^O mod 2}, 
L~ = {uaa + n^b + UcC \ Ua + ni, + Uc — 1 mod 2} + d. 

We say a dimer is part of the n-th layer with n = ria + + ric- The even (odd) layers 
contain positive (negative) dimers. 

The lattice and offset vectors contain 12 coordinate variables, which means the set V is 
a 12-dimensional linear space. However many of the configurations have overlaps, and thus 
are not packings. 

B. Three-parameter family 

In order to verify that the dimers of a double dimer configuration are not overlapping (i.e. 
that it is a packing), intersections of dimers with their neighbors are studied. It is sufficient 
to consider the dimer at the origin: 

Lemma 1. The space group of a double dimer configuration acts transitively on the dimers. 
Proof. The inversion x ^ d + a — x maps positive dimers to negative dimers. □ 
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Intersections can occur between three types of simplices forming the surface of a dimer: 
triangular faces, edges, and vertices. Given two simphces, we say the incidence condition 
for these simphces holds if their intersection is non-empty. Incidence conditions are of either 
face-to-face, face-to-edge, face-to- vertex, edge-to-edge, edge-to-vertex or vertex-to-vertex 
type. 

Next, we replace the simplices by affine subspaces and use the notation V[p] = {p} to 
be the vertex p, E\p,q] = {p + K,{q — p) : n E R} to be the line containing the edge from 
p to q, and F\p, q, r] = {p + K{q — p) + A(r — p) : A e R} to be the plane containing the 
triangular face spanned by p, q, and r. We say the linear incidence condition holds if the 
intersection of two affine subspaces is non-empty. In the following we will consider linear 
incidence conditions only. They are easier to handle analytically and can, in certain cases, 
be used as sufficient conditions for the absence of overlaps. 

Lemma 2. Given the nine linear incidence conditions 

F[o,p,q]n(d-F[o,p,q] + a)^$, 
F[o, q, r]n{d- F[o, q, r]+b)^ 0, 
F[o, r, p]n{d- F[o, r, p] + c) 0, 
F[o,r,p] n {d- F[o,r,p] +a + b + c) 

F[s,r,g]n(ci-F[s,r,g]-o)^0, (3) 
F[s,p,r] n {d- F[s,p,r] -b) j^d), 
F[s,q,p] n (d - F[s,q,p] - c) ^ 
F[s,q,p] Did- F[s, q,p\-a-b-c)^^, 
E[o, q] n {E[s, r] + a + 6) 7^ 0, 

the collection V — {p E V : the conditions G'^,G^,Gf,G^_^_fj_^^,Ga+b hold} is a three- 
dimensional linear space. 

Proof. The nine linear incidence conditions are sufficient for the statement that the dimer 
at the origin is not overlapping with the neighbors at d±a, d±b, d±c, d±{a + b + c), and 
±(a + b). Note that is invariant under inversion at the origin, so that we do not have to 
consider the neighbors at -\-{a-\- b) and — (a + 6) separately. 

Since the planes involved in the face-to-face condition are parallel, they coincide if 
they are incident, so we can alternatively write G^ as an equation: 

G+ : 2o + k{p -o) + X{q-o) =d + a (4) 

with parameters k, A G R. This equation is a vector equation and thus there is one scalar 
equation for each of its three components. Two of the scalar equations are used to eliminate 
K and A. The remaining scalar equation is a single linear constraint for the components 
of the vectors d and a. Similar constraints are obtained for the seven other face-to-face 
conditions. The edge-to-edge condition Ga+b corresponds to the equation 

Ga+b '■ o + ii{q - a) = s ^- v{r - s) + a ^- b (5) 

with parameters G R. Again, a single linear constraint for the components of the 
vectors a and b is obtained. The system of linear constraints turns out to have rank 9, so 
the dimension of the space of double dimer configurations V is reduced from 12 to 3. With 



g: 

G- 
Gb 

g: 

(~<- 

^a+b+c 



5 



the new scalar parameters m, w, the lattice vectors and the offset vector can be written 

as 

^=(-|)-«'i+^'i-2w-^>, 

c = {^-u + 2v + 2w, -f^ + lu - V + 3w,l^ + lu - V + w) , 

□ 

Remark 1. (i) The eight face-to-face incidence conditions are important for dense packings, 
because the corresponding faces are parallel and therefore can be matched perfectly. 

(ii) The face-to-face constraints imposed by G^_^^_,_^ and G~_^j,_,_^ involve the faces spanned 
by o, r, p and s, q, p, respectively. We balance the constraints with an opposing edge-to-edge 
constraint with the neighbors at ±(a -|- b) on the 'opposite' edge spanned by o, q and s, r in 
the form of the incidence condition Ga+b- 



C. Restricted parameter space 

We now specify a subset V" C V of the three-parameter space that consists of double 
dimer packings only, and we later optimize density over this subset. 

Lemma 3. Consider the restricted parameter space 

V" = {{u,v,w) eV : \\u + 2v\ < ^ + w A \v\ < -l^-w} , (7) 

where {u,v,w) is the parameterization in Each double dimer configuration p G V" is a 
packing. 

Proof. Since double dimer packings are transitive, it is sufficient to look for overlaps of 
+F2 only. We write a positive dimer as -I-F2 + n^a + n^b + ricC and a negative dimer as 
d ~ F2 + + Ubb + ricC. The first part of the proof shows that we only need to consider 
the 26 neighbors with —1 < na,niy,nc < 1 when looking for overlaps. Note that V" has 
the shape of a sheared tetrahedron with the four vertices 61^2 = (±3/10, 0, +3/64) and 
63,4 = (=f3/5, ±3/20, —33/320). We verify that for each configuration corresponding to one 
of the four extremal points Cj, the plane a + Hb + Re separates +F2 from +F2 ± 2a and 
the plane d/2 + a + Hb + He separates -I-F2 from d — F2 + 2a. This shows that < 1. By 
similar arguments we get ria > —1 and —1 < nb,nc < 1. Since a, b, c are linear in u, v, 
w and V" is convex, the separating planes and the restrictions —1 < na,nb,nc < 1 remain 
valid for every double dimer configuration in V". 

By construction of V in Lemma [2] we know that the ten dimers c? — F2 ± a, — F2 ± 6, 
d — F2±c, c? — F2±(a + 6 + c), +F2±(ci±&) do not overlap with +F2. Next, we find that the 
dimers +F2 and d — F2 + a + b — c are separated by the plane F[o, p, q] for any configuration 
in V". We find similar planes for d — F2 + a — b + c, d — F2 — a + b + c, d — F2 — a — b + c, 
d — F2 — a + b — c, and d — F2 + a — b — c, which shows that these six dimers cannot overlap 
with +F2. Finally note that pairs of positive neighbors are related by inversion symmetry. 
These considerations reduce the number of neighbors that have to be checked from 26 to 
five. The remaining ones are: +F2 + b + c, +F2 + c + a, +F2 + a — b, +F2 + b — c, +F2 + c — a. 
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Consider the linear edge-to-edge and vertex-to- vertex incidence conditions 



Hc+a 

iib-c 

iir-n 



E[o,r] n {E[s,p]+b + c) ^ 0, 
E[o,p] n {E[s,q]+c + a) ^ 0, 
V[p]r]{E[r,q] + a-b) 
V[q]niF[s,p,r] + b-c) ^ 0, 
V[r] n {F[o,q,p]+c-a) ^ 



(8) 



Similar to the linear incidence conditions in Lemma |2| the conditions tlb+c-, Hc+a, Hf,_c, 
Hc_a can be written as four linear constraints for the parameters m, f , w, which correspond 
to four boundary planes. Each of the planes separates V into two half-spaces. Overlaps can 
only appear in one of the half-spaces, and the half-spaces that are free from overlaps are: 



Hc+a 
H;,_c 

Hc_a 



+ \U + 2V -W < 320, 



V 



\u — 2v 
w < 
+ V + w < 



w < 



33 
320 
33 
320' 



3_ 
64' 
3 



(9) 



The intersection of the four half-spaces is the restricted parameter space V". 

The condition Ha-6 is different from the other ones, since it is a vertex-to-edge condition. 
It is true if and only if both vertex-to-face conditions V[p] fl {F[o,r,q] + a — b) 7^ and 
V\p] n {F[s,r, q] + a — b) (/) are true. We calculate the region of parameter space that is 
free from overlap using the half-spaces for these two conditions: 



H 



a—b 



-M < A +M < 0, 



(10) 



which covers all of V . Hence, the vertex V[p\ can be incident to the edge i?[r, q\ only along 
the central plane m = 0, and tla-b is never needed to prevent overlaps within V' . □ 

Lemma 4. The collection V" contains exactly two maximal density packings with density 
IPj. The packings are related by a crystallographic symmetry operation. 

Proof. The unit cell volume as a function of the parameters {u,v,w), 



V = \det[a + b,b + c,c + a]\ = ^(117 + 60m^ 



80uv - 80v' 



(11) 



is a hyperbolic paraboloid with saddle point at u = v = 0. The parameters {u, v, w) can be 
chosen such that V does not depend on w, which means packings corresponding to the line 
(0, 0, w) are related by a lattice shear. Surfaces of equal volume are hyperbolic cylinders. 
The extrema of V\p" are taken on the boundary. We find that the densest packings are 
located on the maximal lines: 



Hfe+c A Hc_a 
H;,_c A Hc+a 



+ 3|;(2,5,0) + |n(6, 
- 3|(2,5,0) + |m (6, 



-1), 
-!)• 



(12) 



Maximizing the packing density = 2U/V with the tetrahedron volume U = 
^\det[o — p,o — q,o — r]\ = 9 along the maximal lines yields = and two optimal 
points: 



{u±,v±,w±) = i^lo (2,5,0) . 



(13) 
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It can be shown that for two packings with parameters {u,v,w) and {u,v,w) related by 
{u,v,w) = {—u, —v,w) there exists a direct isometry 



T = i -2 -2 +1 (14) 




that symmetrically maps the lattice vectors Ta = —b, Tb = —a, Tc = — c, Td = d and the 
vertices of the dimers To = s, Tp = p, Tq = f, Tr = q, Ts = o. In particular, T maps the 
optimal points onto each other. □ 



D. A dense double dimer packing 

The following theorem is the main result of our work. For the classification of crystallo- 
graphic point groups and space groups we follow the terminology in ^24j . 

Theorem 1. There exists a double dimer packing of tetrahedra with packing density 
Its space group is PI (point group 1 ) and acts transitively on the dimers. 

Proof. The double dimer packing at the optimal point has density It is 

specified by the lattice vectors 

a= Jo (290, 107, -7) , 

^= Jo (-34, 277, 135), (15) 
c= J^(94,-83,247), 

and the offset vector 

d = J^ (38, 5, -25) . (16) 

We know from the angles between the lattice vectors that the crystal system is triclinic. The 
only non-translational symmetries are inversions mapping onto L~, which determines 
the space group. □ 

Remark 2. (i) The lattice spanned by the vectors a, 6, and c in Theorem [l] can be obtained 
from a simple cubic lattice by a small deformation. Furthermore, the norm of d is much 
smaller than the norms of a, 6, c. These two observations mean that the set L = U L~ 
is structurally similar to the rock-salt lattice. 

(ii) One of the axes that is almost a three-fold symmetry is the diagonal a + b + c, which 



coincides roughly with the dimer axis. By the choice of basis vectors in (15) and (16), the 
three-fold symmetry of the individual dimers is broken. 



E. Comparison with previously found double dimer packings 

The one-parameter family of dense tetrahedron packings found by Kallus, Elser and 
Gravel ^22j has constant packing density = ^fy ^ind corresponds to the symmetric line 
(0,0,w;). The packings are transitive on individual tetrahedra as a result of the symmetry 
T, a two-fold rotation around the axis through the vertex p = (2,-1,-1). 

Torquato-Jiao's two-parameter family of packings [2^ corresponds to the generic plane 
bu = —2v. The intersections of this plane with the maximal lines gives the parameters 
2^(± (16, —40, 0) — (0, 0, 5)) and the maximum packing density (p = jflfl- The two packings 
are not related by symmetry. 
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III. NUMERICAL COMPRESSION OF SMALL CELLS 

A. Motivation and methods 

The packing given in Theorem [l] is optimal under two assumptions: (1) the densest 
tetrahedron packing is a double dimer packing, and (2) the packing is a configuration in the 
restricted parameter space V" . To investigate the possibility of denser tetrahedron packings 
- double dimer packings as well as different types of packings - we rely on a numerical 
search. 

We use standard Monte Carlo simulation [25], which allows a system of tetrahedra to find 
dense packings by stochastically exploring all possible configurations subject to the laws of 
statistical mechanics. An elementary simulation step consists of a random displacement 
move within a finite simulation cell - taking into account periodic boundary conditions - 
and a random rotation move of randomly chosen tetrahedra. The move is rejected if it 
generates an overlap, or accepted otherwise. The starting configuration is a dilute, random 
arrangement. During the simulation run, the system is slowly compressed by rescaling the 
size of the simulation cell. The compression is controlled by applying external forces in the 
isobaric ensemble. Additionally, fiuctuations of cell shape by shearing the cell in random 
directions by random amounts are allowed. A lattice reduction technique minimizes the 
distortion of the simulation cell after each shear. For details of the tetrahedron overlap 
detection algorithm we refer to method two in the Methods section of Ref. [?T]. 

The Monte Carlo scheme samples the high-dimensional configuration space stochastically, 
and is not biased towards any particular type of packing. The only constraint is the number 
of tetrahedra in the simulation cell, A^, which does not change during the simulation. In the 
following, the search is restricted to small cells, 1 < < 16, where efficient compressions 
are easily and rapidly possible to high accuracy. 



B. Monte Carlo simulation results 

To obtain sufficient statistics, M = 400 compression simulations were run for each value 
of (1000 runs for = 16). Each run involves 7 million Monte Carlo moves per particle 
and results in a final density (pi, i = 1, . . . ,M for a given A^. The maximum numerical 
density for a given N is cf) = max{0j}. The distribution of (pi close to (p indicates the ease 
with which we can obtain the optimal packing in simulation. For most A^ we find a clear gap 
separating a set of very dense packings from the rest in the sense that the relative density 
variation among the very dense packings is significantly smaller than the gap between very 
dense and less dense packings. In this case, we denote the simulations corresponding to the 
very dense packings as successful. 

The maximum numerical densities are given in Table [T] and their corresponding packings 
are depicted in Figure [1} As can be seen from the success rates in the Table, the geometrically 
constructed optimal packings for A^ = 1, 2 are obtained very easily and in every simulation. 
The A^ = 3 packing is three-fold symmetric and, as far as we know, has not yet been 
reported in the literature. Its optimal density 0=1 can be calculated analytically. The 



structures with A^ = 4, 8, 12, 16 are the dimer packings discussed in section II D The packing 
with N = 5 consists of imperfect pentamers, i.e. four tetrahedra arranged face-to-face to a 
central one. A^ = 6 is a mixture of dimers and single tetrahedra (monomers), and A^ = 7 is 
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#Tetra 


Maximum Density 


Success 


Motifs, 


N 


Numerical, (p 


Analytical, (p 


Rate 


Structural Description 


1 


0.367346 


18/49 


100% 


1 monomer [11] 


2 


0.719486 




100% 


2 monomers, transitive |22j 


3 


0.666665 


2/3 


21% 


3 monomers, three-fold symmetric 


4 


0.856347 


4000/4671 


80% 


2 dimers (positive + negative) 


5 


0.748096 




22% 


1 pentamer, asymmetric 


6 


0.764058 




11% 


2 dimers + 2 monomers 


7 


0.749304 


3500/4671 


15% 


2x2 dimers minus 1 monomer 


8 


0.856347 


4000/4671 


44% 


2x2 dimers, identical to = 4 


9 


0.766081 




— 


1 pentagonal dipyramid + 2 dimers 


10 


0.829282 




2% 


2 pentagonal dipyramids 


11 


0.794604 




— 


1 nonamer + 2 monomers 


12 


0.856347 


4000/ 4o7i 


3% 


3x2 dimers, identical to A = 4 




0.788728 




4% 


1 peiiua^oiiai uipyiaiiiiu. ^ ^ ciimeio 


14 


0.816834 




3% 


2 pentagonal dipyramids + 2 dimers 


15 


0.788693 






Disordered, non-optimal 


16 


0.856342 


4000/4671 


< 1% 


4x2 dimers, identical to A = 4 


8 X 82 


0.850267 






Quasicrystal approximant [21] 



TABLE I: Maximum numerical densities (j) for packings with small cells, obtained with numer- 
ical compression via Monte Carlo compression starting from a random configuration. A data 
file with the packings may be downloaded from the internet at [26]. For comparison, the qua- 
sicrystal approximant result with A^ = 8 x 82 is included. Details about the analytical results 
02 = 9/ (139 - 40\/l0), (p5 = 0.74809657. . ., (pe = 11228544/ (97802181 - 132043^/396129), and 
010 = 29611698560/ (23657426736 + 4919428689^/6) are given in Appendix D of Ref. [23. 

identical to = 8 with one complete vacancy (missing tetrahedron). 

For larger simulation cells crystalline packings are harder to achieve. The = 10 packing 
consists of two pentagonal dipyramids, perfect in the sense that four of each set of five 
tetrahedra are arranged face-to-face; the fifth tetrahedron is oriented in such a way to 
distribute the (obligatory) gap [U] evenly on its two sides. In the cases A^ = 9, 11, and 
15 we found no clear density gap separating very dense packings from less dense packings. 
Either the density gap does not exist, or our simulations were not successful in finding the 
optimal packings. The latter must be the case for A^ = 15, since a dimer packing with a 
vacancy can give a density of = fpf = 0.802825 . . . Instead, the compressions became 
trapped in disordered configurations with a network of pentagonal dipyramids, similar to 
that found in previous simulations [2T]. It is noteworthy that all packings for A^ > 4 achieve 
a density > 0sphere- 
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iV=16 



FIG. 1: (color online) Dense packings with up to 16 tetrahedra per cell as obtained from Monte 
Carlo simulations. The densest dimer packing is observed for = 4, 8, 12, 16. 
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IV. CONCLUSIONS 

We have found the densest known packing of tetrahedra with density = ||^. This result 
was obtained as the optimal solution of a three-parameter family of dimer packings, which is 
a generalization of one-parameter and two-parameter families of packings recently reported 
with lower maximum densities. Interestingly, the densest packing is not the most symmetric 
one. Isobaric Monte Carlo simulations with variable cell shape starting from random initial 
conditions recover the same high packing density within 10~^ for small systems containing 

4, 8, 12, and 16 tetrahedra. We also discovered new candidates for densest packings with 3, 

5, 6, and 10 tetrahedra per unit cell. The analytical and numerical results combined suggest 
that the packing density reported here could be the highest achievable, at least for small A^. 

The dimer structures are remarkable in the relative simplicity of the 4-tetrahedron unit 
cell as compared to the 82-tetrahedron unit cell of the quasicrystal approximant [H] , whose 
density is only slightly less than that of the densest dimer packing. The dodecagonal qua- 
sicrystal is the only ordered phase observed to form from random initial configurations of 
large collections of tetrahedra at moderate densities. It is thus interesting to note that for 
some certain values of N, when the small systems do not form the dimer lattice packing, 
they instead prefer clusters (motifs) present in the quasicrystal and its approximant, pre- 
dominantly pentagonal dipyramids. This suggests that the two types of packings - the dimer 
crystal and the quasicrystal/approximant - may compete, raising interesting questions about 
the relative stability of the two very different structures at finite pressure. These questions 
will be explored in a forthcoming paper. 
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Appendix A: Pictures and equations 

This appendix contains figures of tlie double dimer packings (Figures A1-A3), tlie tliree- 
dimensional restricted parameter space (Figure A4), and tfie volume contour function (figure 
A5), as well as visualizations of the internal symmetry of the restricted parameter space 
(Figures A6 and A7). 

Figure AO introduces the orientation of the axes for all following figures. Wc distinguish 
between the vector coordinate axes (x.y.z) and the packing coordinate axes {u,v,w). In 
the former case, the viewpoints are (1, 1, 1), (O, — cos j^, sin j^), (0, 0, 1), and the ranges are 
< x,y,z < +y (Figures A1-A3) or < x,y,z < +^ (Figure A6). In the latter 
case, the viewpoints are (cos || cos f^, cos || sin , sin || ) , (sin j^, 0, cos and the range 
is -I <u,v,w < +1 (Figures A4-A5 and A7). 

Figure Al (4 pages) gives the notation and basic equations, and shows the nearest inter- 
secting neighbors. The colors of the equations match the color of the respective neighbors. 
Since the double dimer packings are cluster transitive on dimers, we consider only the dimer 
at the origin. Its positive and negative neighbors are shown for three pairs of special pack- 
ings at the optimal points, the central points, and the symmetric points in parameter space 
(Figure A4). The pairs of packings are related by the symmetry T (Figure A7). Notice that 
the sets of negative neighbors are always the same, whereas the sets of positive neighbors 
are different. The neighbors for the optimal packings are related by a reflection through the 
plane +2x — y — z — 0, the neighbors for the central packings by a reflection through the 
plane +2x — y — z — 0, and the neighbors for the symmetric packings by a reflection through 
the plane which contains ±a ± b. 

Figure A2 (2 pages) lists all 20 possible neighbors that can intersect with the dimer at 
the origin together with their linear incidence conditions. The entire figure is a floor plan 
diagram of the neighbors, ordered by layer. Notice how the positive layers and the negative 
layers complement each other. The even layers are related by symmetry, which is broken for 
the odd layers due to the presence of a small translation by the vector d. 

Figure A3 (2 pages) shows larger portions of the layers for the optimal packings. The 
entire flgure is a floor plan diagram of the packing, ordered by layer. Notice the almost 
triangular/hexagonal planar basis mentioned in Remarks 2. We show both optimal packings 
for comparison to demonstrate the subtle differences. The symmetry T that maps the 
packings onto each other is a two-fold rotation around an axis running from bottom left to 
top right in the left columns (view along the 'almost' three-fold axis). 

Figure A4 (1 page) depicts the restricted three-dimensional parameter space of dense dou- 
ble dimer packings {u, v, w) together with special planes, lines, and points. The shape is a 
strongly deformed tetrahedron extended along the -u-axis. It is bounded by four boundary 
planes. Packings with maximal (minimal) density are located on the maximal (minimal) 
lines (Figure A5). The optimal plane contains the two optimal points, the origin, and the 
symmetric line. 

Figure A5 (1 page) visuahzes contours of the unit cell volume V — ^{in + dOu^ — SOuv — 
80v^) in parameter space. The lattice volume is a hyperbolic paraboloid in {u,v,V). The 
contours are hyperbolic cylinders in {u,v,w) parallel to the w-axis. Their intersection with 
the restricted parameter space (gray) allows to identify the locations of volume extrema on 
the maximal and minimal lines. 
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Figure A6 (2 pages) demonstrates the distortion in the geometry of nearest neighbors in 
vector space {x,y,z). The colored spheres (colors match the colors in Figures A1-A3) are 
positioned at the centers of the twelve dimers in the layer 0, ±1, closest to the dimer at 
the origin. Lines connect dimers that are adjacent in the dimer lattice. A deformation of 
the neighbor shell occurs when moving along the various lines in parameter space {u,v,w). 
This deformation is illustrated in each sub-figure by overlaying three neighbor configurations 
along these lines. Packings along the central line have the symmetry T, which is a two-fold 
rotation about the axis p = (-1-2, —1, —1) (see proof of Theorem 2). Notice how these 
packings and distortions look more symmetric than the others. Packings along the maximal 
lines form a double family. Packings along the minimal lines form two single families that 
are distorted in complementary ways. 

Figure A7 (1 page) illustrates the effect of the symmetry T on the restricted parameter 
space shown in Figure A4 by identifying the equivalent packings under this symmetry: 
T {u,v,w) = {—u, —v,w). Within contour planes of the form w = *, i.e. normal to the w- 
axis, T acts as a combination of a shear and a two-fold rotation around the symmetric line 
(0, 0, w). Equivalent points are colored identically. Along the minimal lines, opposite points 
on the same line are related. Notice that each minimal line is contained in a single contour 
plane. Along the maximal lines, opposite points on opposite lines are related. Notice that 
the two optimal points are contained in a single contour plane. The final sub-figure shows 
the change of the contour plane intersection with the maximal lines for a sequence of five 
contour planes. 




Fig. AOvec The vector {x,y,z.) coordinate axes for Figs. Al - A3 & A6 




Fig. AOpac The packing {u,v,w) coordinate axes for Figs. A4-A5 & A7 
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The cluster F2 has vertices 

o = (+2,+2,+2) 

p = (+2-\-\) 

q = {-\,+2-\) 

r = {-\-\,+2) 

s = (-2-2-2) 
The cluster volume is 

U = ^det[o-p,o-q,o-r]+^d&t[s-p,s-q,s-r] = 18 
Clusters have 2 possible orientations 

+F2 positive 

-F2 negative (point reflection = scalar multiplication by - 1 ) 
The family of packings is parametized by 

27 51 753 1 

{a^,by,c,} = <+-,+-,+ -y)+{H,v,-M-v+w) 
The lattice vectors are 

a = {a^,ay,a.) = <+^,+|-|,>+h< 1,0,2)- v<0,+ 1,-1) 
b = (b,,b,A) = <-fo,+|'+2-o>-«<l'0'2)+v<0,+ l,-l) 
c = {c„Cy,c,) = {+\^-^£i,+]fQ}-(jU-v)(+2-l-l}+w{2,3,l} 

d = (d,,dyA> = ( + Po,-2-o,-4)+«<l'l'l>+^<0,+ l,-l> 

The lattice volume and packing density are 
V = 2det[a,ii,<:] = det[a+b,b+c,c+a] 



^ 2U/V - 36/V 




Fig. Al|^ , The 14 intersecting neighbors of the optimal packing {u,v,w) = {+ T,,+i-^0) 
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The cluster +F2 has 20 neighbors in various layers 
The 1 2 positive neighbors are 

+¥2+a+h +F2+b+c 



+F2+a-h 
+¥2-a+h 
+F2-a-b 



+F2+/5-C 
+F2-b+c 
+F2-b-c 



The 8 negative neighbors are 

-F2+d+a+b+c 



-F2+d+a 



-F2+cl+b 
-Fo+d-b 
-F2+d-a-b-c 



+F2 +c+a 
+F2+c-a 
+F2-c+a 
^Fj-c-a 



~F2+d+c 

-F2+d-c 



The 2 optimal packings are parametrized by 



The lattice vectors are 
21 



a 



, 87 3 2_l _ ^ , 



^~160'"'"320'"'"64^ 
_ , 141 249 741, 

^ ~ ^"'"l60'~32b'"'"320^ 



^"""160' 64' 



The lattice volume and packing density are 

_ 42039 

~ 1000 
4000 

^ ^ 4-671 



+2 layer 
layer 
layer 

-2 layer 

+3 layer 
+ 1 layer 
-1 layer 
-3 layer 



^ " ^"^160' "'"320' 64' 

h-(-- ^- 

~ '~32'"'"3l0'"'"32b^ 

^"^"'"160' 64'"*" 64^ 
, _ , 13 _37 _ 7 , 

^^"""160' 320' 320^ 



. . _ 42039 

~ 1000 
4000 

^ ^ 4671 




Fig. Al„ ( The 14 intersecting neighbors of the optimal packing {u,v,w) = {- 



160' 64' 
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g: 



g: 





































/ 


\ 



H„_„ V{p^r\E{r,q\+a-b 
Hft_, V[ci\r\F[s,p,r]+b-c 
He -a V[r] n F[o,q,p\+c-a 



Yi n F2+a-b +x-y-z = 4 
Ft n F2+b-c 



+5x-y-z = 12 
-x+5y-z = 12 



F2 n F2+c-a 



-y+z = 4 




F[o,g,r] n d-F[o,q,r]+b 
F[o,r,p] n £/-f [o,r,/7]+c 



F2r\d-F2+a 
F2 n d-F2+b 
Fi n ^/-^2+c 



Hft+, E[o,r]r\E[s,p]+b+c 
He +a E[o,p] n E[s,q\ +c+a 



Fj n F2+fl+Z> 

F2 n Fj+^+c 

F2 n F2+c+a 



+x+y-z - 2 
-.v+>'+z = 2 
+x-_y+2 = 2 




+x+3y-z = 6 
-x+>'+3z = 6 
+3x-v+z = 6 




+3 layer 

Ga+i,+e-f[o,'",p] n ii-F[o,r,/7]+fl+Z7+c F2 n d-Fj+a+b+c +x-y+z = 2 
Fig. A2^ The neighbors & intersection equations & separating planes of F2 (positive layers) 
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H„_„ E[r,q\r\V[p]-a+b 
n„-c F[s,pA n V[q\-b+c 
^ a F{o,q,p\r\V[r]-c+a 



¥2 n ¥2-a+b -x+y+z = 2 -5x+y+z = 6 
F2 n F2-b+c +x-5y+z = 6 

F2 n Fn-c+a +x+y-z = 2 




F[s,r,q]r\d-F[s,r,q]-a F2r\d-F2-a 
F[s,p,r]nd-F[s,p,r]-b FzH^j'-rz-^ 
GJ F[s,q,p]nd-F[s,q,p]-c F^flfZ-F.-c 




Ga+b E[s,r] n E[o,q]-a-b 
Ub+c E[s,p] n E[o,r]-b-c 
Uc +a E[s,q] n E[o,p] -c-a 



¥2 n ¥2-a-b 
¥2 n F2-/?-c 
FzHFz-c-a 




-5x+y+z - 6 
+x-5)'+z = 6 
+x+y-5z = 6 





£i 











-x-3y+z = 6 
+x-)'-3z = 6 
-3x+v-z = 6 




-3 layer 

Q~i^j^^,^^F\s,q,p\ n d-F\s,q,p\—a-b-c ¥2 fl d-¥2-a-b-c +x+y-5z = 6 
Fig. A2~ The neighbors & intersection equations & separating planes of F2 (negative layers) 
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optimal plane 
central plane 
boundary plane 
boundary plane 
boundary plane 
boundary plane 



symmetric line 
minimal line 
minimal line 
maximal line 
maximal line 



origin <0,0,0> 
optimal point 
optimal point 
central point 
central point 
symmetric point 
symmetric point 



dV 




H,_, A dV 
Hb-c A H(._„ 

Hf,_(. A Hc+a 



H/,H-. A H,_„ A dV 
Hi_, A H,+„ A dV 

Hh+r A Hf-n A H„ 



H/,-r A Hc+a 

B-h-c A H(._(, 
H,,, A H 



1 ^ 33 
I ^ 33 

--U-2V-W = 3- 

3 

-v+w = - 

6,4 

+ V+VF = - 


















> / 


Y / \\ 



(0,0,0)+R'(0,0,l> 



<0,0-3|>+M<l -^,0> 




A H„_,, A (5V 
A H„_fr A 5V 



(- 0) 

^ 160^' 64'"' 

(0 -1 - 1 > 

^ ' 20' 320' 

(0,0,-|„> 



Fig. A4 Special planes & lines & points of the parameter space {u,v,w) 



V = |(117+60i<--80Hy-80v-) 



<p = 100/(1 17+60M^-80wy-80v-, 



^ _ 42093 

1000 
, 4000 

^ = 4-671 

maximal lines 

tangent to contour 
optimal plane 

contains tangent points 
optimal points 
= tangent points 



1053 42093 

— > V > — 

25 1000 
100 , 401X1 

Tl7 < < 4-671 



V = 1^'-' 

25 

, 100 

'^"Til 

asymptotes of contour 

= contour 
symmetric line 

= central axis of contour 
symmetric points 

@ central axis of contour 


















V 












Fig. A5 Intersections of the parameter space {u,v,w) with contours of lattice volume V 



^=1 



+ 1 


-2 


-2 


-2 


-2 


+ 1 


-2 


+ 1 


-2 



T : d,a,b,c -> d,-b,-a,-c 
fixed line /? = (+2,-1,-1) 
in vector space (jc,y,z) 

any 2 packings related by 

(¥,V,M7) — V,+w) 

fixed line (0,0,0) +w(0,0,l). 
fixed planes w = * 
has symmetry T 

symmetric line 
sym to self 
origin 

@ w = 

sym to self 

minimal point 

3 

@ w - + — 
64 

sym to self 
minimal point 

@ VV = 

320 

sym to self 



3 

+6-4 




minimal line 

@ w ■ 

sym to opposite point 
@ same minimal line 
minimal line 

^ 33 

@ VV = 

320 

sym to opposite point 
@ same minimal line 

2 maximal lines 

skew lines, no single plane 

spans family of planes 
sym to opposite point 

@ opposite maximal line 
2 optimal points 

@ w = 

sym to each other 
2 central points 
@ w = — — 

^ 320 

sym to each other 
Fig. A7 Packings having symmetry 7", or packing pairs related by the symmetry T 
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Appendix B: Symmetry 
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BO Origin (0, 0, 0) 

The symmetry group is Z| 



10 
10 
1 



-'o 



Tq : packing (0, 0, 0) ^ packing (0, 0, 0) 
d,a,b,c d,+a,+b,+c 
o,p,q,r,s o,p,q,r,s 



c 624 



+538 +172 -344 
+79 +466 +316 
-251 +502 -380 



: packing (0, 0, 0) —?■ packing (0, 0, 0) 
d, a, b, c d, +a, +b, —c 



ab 208 



+98 -196 -24 
-165 -86 -36 
-55 -98 +196 



rp-l 



T„(, : packing (0, 0, 0) -> packing (0, 0, 0) 
d, a, b, c d, —a, —b, +c 



T = i 

abc 3 



+1 -2 -2 

-2 -2 +1 
-2 +1 -2 



^abc rotation by ±7r about p — (+2, —1,-1) 



r„(,^ : packing (0,0,0) packing (0,0,0) 
d, a, 6, c d, —a, —6, — c 

o,p,q,r,s s,p,r,q,o 

Bl Symmetric line (0, 0, 0) + w (0, 0, 1) 

The symmetry group is Z 



S 



iw) 



o-l 



10 




+2 -4 +8 


10 


+3 -6 +12 


1 


_+l -2 +4_ 


packing (0, 0, 0) — )■ pack] 


d, a 


,b,c 


d, a, b, 


10 


117 


+2 -4 +8 


10 


+3 -6 +12 


1 


_+l -2 +4_ 



rr-l 



5, 



S^^^ : packing (0, 0, 0) packing (0, 0, -w) 



d, a, b, c 



d, a, b, c 
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B2 Symmetric line (0, 0,0) +w (0, 0, 1) 
The symmetry group is Z| 

^(w)'^oSi'w) ■ packing (0,0, w) packing (0,0, w) 



d, a, b, c 
o,p,q,r,s 



d, +a, +b, +c 
o,p,q,r,s 



20 
117 



W 



+2 -4 +8 

+3 -6 +12 
+1 -2 +4 



^{w)TcS(^) ■ packing (0,0, w) -> packing (0,0, w) 



d, a, 6, c 



d, +a, +6, — c 



20 

ab 117 



+2 -4 +8 

+3 -6 +12 
+1 -2 +4 



-^H^afc-^H • packing (0, 0, w) packing (0, 0, w) 
d,a,b,c d,—a,—b,+c 

S(w)'^abcS{2) = ^abc rotation by ±7r about p = (+2, -1, -1) 
^(w)'^abcS{w) '■ packing (0, 0, w) packing (0, 0, w) 



d, —a, —b, —c 
s,p, r,q,o 



d, a, b, c 
o,p,q,r,s 

B3 Pcirameter space 

For any packing {u, v, w) 

— : {x, y,z) — > d — {x, y,z) + a inversion about |(d + a) 

i?„ : packing (u, v, w) — >■ packing {u, v, w) 
d,a,b,c a,d,d + a — b,d — c + a 

Rf^ — R^^ : {x,y,z) — >■ d—{x,y,z)+b inversion about |(d + 6) 

Rf^ : packing {u,v,w) packing (u.v.w) 
d,b,c,a b,d,d + b — c, d — a + b 

R^ — R~^ : {x, y,z) — >■ d — (x, y,z) + c inversion about |(d + c) 

R^ : packing {u,v,w) — >■ packing (u.v.w) 
d,c,a,b c,d,d + c — a,d — b + c 

For any 2 packings {u,v,w), {u,v,w) related by 

(■u, v,—^u + V + w) — {—u, —V, —\u^-v^-w) 



T = i 

ahc 3 



+1 -2 -2 

-2 -2 +1 
-2 +1 -2 



rotation by ±7r about p = (+2, —1, —1) 



^abc • packing {u,v,w) — )■ packing {u,v,w) 
d,a,b,c d,—a,—b,—c 
o,p,q,r,s s,p,r,q,o 
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Appendix C: Computations 

CO^gp Vector space {x, y, z) 

The cluster vertices are 

0= (+2, +2, +2) 
p= (+2,-1,-1) 
g = (-l,+2,-l) 
r = (-l,-l,+2) 
s = (-2,-2,-2) 

The cluster volume is 

[/ = I det[o — p, o — g, o — r] + I det[s — p, s — g, s — r] = 18 

The lattice vectors are 

The lattice volume and packing density are 

V = 2 det[a, b, c] = det[a + b,b + c,c + a] 
= 2U/V = 36/y 

COpj^c Packing parameter space {u, v, w) 

The linear incidence conditions are 
r^± /~(± /-(it /-(± /-( 

The lattice vectors are 

a = {%+u,§-v,-^ + 2u + v) 

c = - 2x4 + 4^; + 2w, -|f - M + 2^; + 3w, g§ + 

^ = ( ]^ + ^' + ^ + ^' + ^ - ^> 
The lattice volume and packing density are 

F = 2det[a,6,c] = |(117 + GOii^ - gOw^; - 80^;^) 
= 36/y = 100/(117 + 60x^2 - 80uv - 80v^) 

Clgy/Clgp^ Optimal plane u—^v 

The lattice volume function is a rotated parabolic hypcrboloid 

V = 1(117 + 60^2 - 80uv - 80v^) 
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The contours (of constant V) are rotated hyperbolas 

5 5 5 36 100 

The asymptotes (of the hyperbolas) are 
u = +2v 

The 2 maximal lines have slope du/dv — —6 

/ + ^ I A _A\ I ^/i _1 I 5\ 

\^ 160 '^64' 80/ ^ \ ' 6' ^6/ 
\ 160' 64'~80/~"'\ ' 6'~2/ 

The contours have slope du/dv — —6 when 
u = +lv 

Cl„_{,/Clj,(,jj Central plane m = 
The linear incidence conditions are 

G± /-^± /-^± r^± ( ' I I 
a ' ' ' *^abc' "^a+fe' "tla-b 

The lattice vectors are 

a = /27 21 _ _A +^;\ 
" \10' 20 '^' 20 ^ ^/ 
5 = /_A 51 ^ 27 _ \ 
^ \ 10' 20 ~ '^' 20 ^ / 

c=(ii+4t' + 2«^'-l + 2t' + 3«;,l+«;> 

(t0'^20"'~'^'~20~'^) 

The lattice volume and packing density are 

y = ^(117-80^2) 
= 100/(117 - 80^;2) 

Cl^_,_^ Boundary plane + v — w — ^ 

The linear incidence conditions are 

(~\± (~i± (~i I I 

' ' ' ^abC ^a+fe' -"-b+c 

The lattice vectors are 

« = (i + ^'i-^'-|) + 2M + ^> 
^ = <-i!)-«'i + ^'i-2«-^^> 

c = (I + 6v, + 2u + 5v, I + M + v) 

d= (Y^ + u,-j^+u + v,-i^ + u-v) 

The lattice volume and packing density are 

V = ^(117 + 60^2 - 8Quv - 80^2) 
= 100/(117 + 60x^2 - 80uv - 80v'^) 
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Clg_j.„ Boundeiry plane —3v — w — ^ 

The linear incidence conditions are 

' G^b^' G^afcc' H^+a 

The lattice vectors are 

c= {l-2u-2v,-%-u-7v,l- 3v) 
d ^ {jq+u,—^+u + v,—^ + u-v) 

The lattice volume and packing density are 

V = 1(117 + 60^2 - 80uv - SOv^) 

(j) = 100/(117 + 60^2 - 80mv - 80^2) 

Clj_^ Boundeiry plane —u -\-2w — ^ 
The linear incidence conditions are 

/-(it /~<± /-(± /-~(± /-( TT 

i^h i^c ^ ^abci ^a+bi "-b-c 

{a,,hy)^{u+%,v + %) 
The lattice vectors are 

c = ( A - « + 4t;, -I + + 2t;, f + 
(i = ( + + M + u, + M - 

The lattice volume and packing density are 

V = 1(117 + 60^2 - SOuv - 80^2) 

= 100/(117 + 60^2 - SOuv - 80v^) 

Clg_„ Bound£try plane —u + Av + 2w — ^ 
The linear incidence conditions are 

) G^, G^, G^f^^, G^_,_^, H^_^ 

The lattice vectors are 

« = <i + «'i-^'-|) + 2^^ + ^^> 

= + M, + U + V,-^+U- v) 

The lattice volume and packing density are 

V = 1(117 + 60?/^ - 80Mf - 80f2) 

(j) = 100/(117 + 60^2 - SOuv - 80v^) 
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C2,y^ Symmetric line (0, 0,0) +w (0, 0, 1) 

The linear incidence conditions are 

, G^, G^, G^^^, G„^(„ H^_6, dV 
{by,c,)^{v+%,w+'^),v^O[mm V] 

The lattice vectors are 

a = /27 21 _3_\ 
" \ 10' 20' 20/ 
b= /-A 51 27\ 
" \ 10' 20' 20/ 

\lO' 20' 20/ 

The lattice volume and packing density are 

T/ _ 1053 
^ 25 

^ 117 

C2+i^ Minimal line (0, 0, + u (l, 0, +|> 
The linear incidence conditions are 

Ga , G^, G^, G^^^, G„_^j„ Hj_^, H^_„ 

The lattice vectors are 

« = <i + «'i'-|) + 2«> 

" \ 10 "'20' 20 ^"/ 

c = (l5-«'-i + ^«'f + i«> 

^= (]^+^'-^+^'-2b+^> 

The lattice volume and packing density are 

V = 1(117 + 60^2) 
(j) = 100/(117 + 60^2) 

02-;^ Minimal line {0,0,-^) + u (l, +|> 

The linear incidence conditions are 

Ga ) G^, G^, G^j,^, Ga+b^ -^b+C "^c+a 
«x = + 10 

The lattice vectors are 



20 ^ 4 " ' 



« = <i + «'i + i«'- 

" \ 10 "' 20 4"' 20 4"/ 



^ \5 2"' 20 



rf = (^ + M,-^ + |7i,-|j + |n) 
The lattice volume and packing density are 

= 1(117 + 75^2) 
= 100/(117 + 75^2) 
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C2+^ Maximal line (O, +1^,-^)+'^ "I' +§> 
The linear incidence conditions are 



G}^, G^, G^^^, G„_^{„ Hj_^^, Hc-a 

The lattice vectors are 

" \ 10 "'5 e"' 10 6 "/ 
^ \10 5 ^ 6"' 10 ^ 6"/ 

ci = (i + „,0+|«,-^ + |«> 
The lattice volume and packing density are 

^=Hlo(4671 + ^(«-Tion 
= 4000/(4671 + ^{u - 4)2) 

Maximal line {^.-yo^+Ii) + u{l,-\,+\) 
The linear incidence conditions are 

G}^, G;^, G^j, , G„^j, Hj_c, H^+a 

27 

= + 10 
The lattice vectors are 

" \ 10 "'2 e"' 5 6 "/ 
C = /X _ 5 X , 1 12 , 1 \ 

\10 3"' 10 ^ 6"' 5 ^ 2"/ 

The lattice volume and packing density are 

^=Hlo(4671 + ^(.+ Tio)^) 
= 4000/(4671 + ^(« + 4)2) 

C3„,ig Origin (0,0,0) 

The linear incidence conditions are 

G^feC ^a+h 

(a^, 6^, c^> = (li + f^, V + |, ^/; + 11) , (w, v, w) = (0, 0, 0) 

The lattice vectors are 
a = /+27 ,21 _ 3_\ 

" \~10'~20' 20/ 
6= /-A +51 ^27\ 

\ 10' '^20' '^20/ 
r= /-4-i29 _237 I 753 \ 

\~^160' 320'~'~320/ 

\^10' 20' 20/ 

The lattice volume and packing density are 

1053 
25 



V = 2det[a,6,c] - 



36/1^ = 
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C3+ 1 Optimal point (+^, -|> 
The linear incidence conditions are 

G}^, G^, G^^^, G„_^{„ H(,_^^, Hp_a, dV 

The lattice vectors are 



-ib> = (+2.71875, +1.003125, -0.065625) 
^= (-i)'+ii'+i> = (-0.31875, +2.596875, +1.265625) 



c 



141 _ 249 I 741 
160 ' 320 ' 



+ = (+0.88125, -0.778125, +2.315625) 
^ = (+ii' = (+0.11875, +0.015625, -0.078125) 

The lattice volume and packing density are 

V = ^ = 42.039 

= ~ .856347677156 

C3opt Optimal point (- ifo, + ^> 
The linear incidence conditions are 



G^a ' GJ^ , , G^^^, G„^j, H^_c, H^+a) 
= ""+ -lio [n^in ^] 

The lattice vectors are 

« = +11' -i> = (+2.68125, +1.096875, -0.234375) 

6= (-|,+|1,+||) = (-0.28125, +2.503125, +1.434375) 

c = <+Ti' -i' +lf > = (+0.73125, -0.703125, +2.390625) 
d = l+j^, = (+0.08125, -0.115625, -0.021875) 

The lattice volume and packing density are 

V = ^ = 42.039 

= ~ .856347677156 

C3+, Central point (O, +^, -i^) 
The linear incidence conditions are 

' GJ^, G^, G^bc^ "^b+C ^c-a^ "^a-b 

The lattice vectors are 

« = <+i'+i'-n)> 

6= (-A +13 +13\ 
\ 10 ' ^ 5 ' ^ 10 / 

\^10' 5' ^10/ 

^=<+B'0,-^> 
The lattice volume and packing density are 

T/ _ 5256 
^ 125 







125 
146 



Central point (O, - ^ , +^ > 

The linear incidence conditions are 

) G^, G^, G^^^, G„^(„ H^_^, H^+„, H„_ 

The lattice vectors are 

a = /+27 ,11 _i\ 

" X^IO'^IO' 5/ 
^ \ 10'^2'^5/ 
^ \^10' 10' ^ 5 / 

The lattice volume and packing density are 

Y — 5256 
^ 125 

^ 146 

C3+^ Symmetric point (0,0,+^) 
The linear incidence conditions are 

) G^, G^, G^f,^, G„_,_(,, H(,_^, H^_a, 
The lattice vectors are 

a = /+27 ,21 _3_\ 
\^10'^20' 20/ 
5= /_A +51 ^27\ 
" \ 10' ^20' ^20/ 

^ \~10' 5' ~ 5 / 

" \~10' 20' 20/ 

The lattice volume and packing density are 

T/ _ 1053 
^ 25 

^ 117 

C3-y^ Symmetric point (0,0, 

The linear incidence conditions are 

± ± (~\± f ' M M M 

) ? ) ^abc^ ^a+b^ "-b+c^ "-c+a^ "^a- 

The lattice vectors are 



a = 


/+27 
\^ 10' 


+ 21 _ 
^20' 


3 

'20 


b = 


(-^ 


+ 51 + 
~20' ^ 


27' 




\ 10' 


20, 


c = 


<+!'- 


_21 +! 
20' 


I) 


d = 


(+1^' 


1 

20' 


1 

'20 



The lattice volume and packing density are 



1053 



V 25 
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Appendix D: Small unit cell packings 



Dl The clusters are 

+B^^+B[o,p,q,r] 
The cluster vertices are 

o=(+l,+l,+l) 

p= (+1.-1.-1) 
g = (-l,+l,-l) 
r = (-l,-l,+l) 

The neighbors of arc 

B^ + a Bi + 6 B^ + c 

B^-a B^-b B^-c 

B^ + a + b B^ + b + c B^ + c + a 

B^ — a — b B^ — b — c B^ — c — a 

B^ + a + b + c 
B^ — a — b — c 

The linear incidence conditions are 

E[o,p] n {E[q,r] + a) ^ 
E[o,g] n (E[r,p] +6) ^ 
E[o,r] n {E[p,q] + c) ^ 
F[o,p,q]n{V[r] + a + b)^0 
F[o,q,r] n {V[p]+b + c) ^ 
F[o,r,p] n {V[q] + c + a) ^ 
V[o] n {F[p, g, r] + a + 6 + c) 7^ 

optimize (3 free parameters) 

The lattice vectors are 



The lattice volume and packing density are 



« = (+2>-|>-i> 

c=<-i'-|'+2> 




196 
27 



18 
49 



D2 The clusters are 



+Bi = +B[o,p,q,r] 
-Bi = -B[o,p,q,r] 



The cluster vertices are 



0= (+1.+1.+1) 
p= (+1,-1,-1) 
g=(-l,+l,-l) 
r = (-l,-l,+l) 



The neighbors of are 

B^ + a + b B^ + b + c B^ + c + a 

B^ — a — b B^ — b — c B^ — c — a 

B^ + 2a + b + c B^ + 2b + c + a Bi + 2c + a + b 

Bi-2a-b-c Bi-26-c-a B^ - 2c - a - b 

Bi + 2a + 2b + 2c 

Bi - 2a - 26 - 2c 

d-B^ + a d-B^ + b d-B^ + c 
d — B^ — a d — B^ — b d — B^ — c 
d-B^ + a + b + c 
d — B^ — a — b — c 

The hnear incidence conditions are 

E[o,r] n {E[p,q] + a + b) ^ 
E[o,p] n {E[q,r] + b + c) 7^ 
E[o,q]n{E[r,p]+c + a)y^(l) 
F[o, q, r] n {V[p] + 2a + 6 + c) ^ 
E[o, r, p] n{V[q]+2b + c + a) 
F[o,p, q] n lv[r] + 2c + a + 6) 7^ 
V[o] n {F\p, q, r] + 2a + 26 + 2c) ^ 

F[o, q, r]n{d- F[o, g, r] + a) 7^ 
F[o,p,q]n{d-F[o,p,q]-b)^(l> 
F[o, r, p]n{d- F[o, r, p] - c) 7^ 
F[o,p,q] n {d-F[o,p,q] + a + 6 + c)^0 
F\p, q, r]n {d — F[p, q,r] — a — b — c) 

optimize (2 free parameters) 

The lattice vectors arc 

a = l (-7 + -19 + 7^, +16 - 4^) 
b=l{+7- iVlO, -7 + iVlO, -10 + A^/lO) 
c = I (-1 + Iv^, +25 - 7V1O, +2 - 2V1O) 
d=l {-l + IVW, +9 - 3v^, -8 + 2V10) 

The lattice volume and packing density are 

= 2det[a,6,c] = det[a + 6, 6 + c, c + a] = p(139-40 
(f) = 2-yv = 9/(139 - 40V^) 

The clusters are 

+Bi = +B[o,p,q,r] 
-Bi = -B[o,p,q,r] 
TBi = TB[o,p,q,r] 
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The cluster vertices are 



0= 

p= (+1.-1.-1) 
q = {-!,+!, -1) 
r = (-l,-l,+l) 



3 



+1 -2 -2 
-2 +1 -2 
-2 -2 +1 



reflection through the plane +x + y + z — 



The neighbors of B^^ are 

B^ + a-b B^ + b-c 



Bi + c — a 
B^-b + c B^-c + a 



a 

B^~a + b 
B^ + 2d 
Bi-2d 

-Bi - a 

TB^ 
TBi 

The linear incidence conditions are 



d + a TB^ + d + b TB^ + d + c 
d + a TBi -d + b TB^ -d + c 



E[r,p]n{V[q\+a-b)^^ 
E\p,q]n{V[r] + b-c)y^(D 
E[q,r]n{V[p] + c-a) ^0 
V[o]n{F[p,q,r]+2d) ^0 

F[o, q, r] n {-F[o, g, r] - a) 7^ 
F[o,r,p] n {-F[o,r,p] - 6) 7^ 
F[o,p,q]r\{-F[o,p,q]-c)^(/} 

E[o,p] n iTE[o,q]+d + a) 7^ 
E[o,q]n{TE[o,r]+d + b)y^tD 
E[o,r] n (TE[o,p] +d + c) ^ 
F[p, q, r] n iTF[p, g, r] - + a) ^ 
F[p, g, r] n ItF\p, q,r] - d + b) ^ $ 
F[p, q, r] n {TF[p, q,r]-d + c) 

The lattice vectors are 

a = (+1,-1,0) 
6 -(0, +1,-1) 
c=(-l,0,+l) 

(+1+1+1) 

The lattice volume and packing density are 

y = 2 det[d + a, d + 6, d + c] = 12 







/y 



The clusters are 

+B5 = +(F2UBpUB^UB,) 

F^^ B[o,p,q,r]UB[so,p,q,r] 
Bp ^ d + TpB[sp,q,r,o] 
Bg = e + T^B[sq,r,o,p] 
B^ = f + T^B[sr,o,p, q] 

The cluster vertices are 



o 
p 

q 

r 



(+1,-1,-1) 
(-1,-1, +1) 



s — 



rp0 

^ {i,j,k) 



(cos( 



1 








+ (sin e) 


-fc -hj 


+ (1 -cos^) 


ii 


ji 


ki 





1 





+k -i 


ij 
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kj 








1 




-j +i 0_ 




ik 


jk 


kk 



T 

Q 



rotation by 9 about the unit vector (i, j, k) 

(-1,+1,+1>/V3 

_ 

(+1 -l,+l>/^/3 
" ^(+l,+l-l>/y3 



The neighbors of Bg are 
B. + a B. + b B 



B. 



B, -6 B 



B^ + a-b Bg 
B.-a + b B. 



+ c 

+ b-C Bg 

-b + c B. 



+ c — a 
— c + a 



The linear incidence conditions are 

F[o,q,r]nid + TpF[o,q,r])^0 
F[o,r,p]n(e + T^F[o,r,p])^0 
F[o,p,q]n{f + T^F[o,p,q])^<I} 

{f + T,E[o,q])n{d + T^E[o,q])j^(!) 
{d + TpE[q,r])n{e + T^E[o,r])^(I} 
{e + T^E[r,o])n{f + T^E[o,p])^iJ) 

E[o,p] n {E[so,q]+a) ^ 

{d + TpE[sp, q]) n {E[so, r] + b)^0 

{e + T^F[sq,o,p])n{V[so] + a)^$ 

{f + T,E[o,q])n{E[so,p]+b)^(Ii 

{d + TpE[sp, r]) n {E[so, g] + c) ^ 

(e + T^E[sq, o\) n (/ + T,E[sr,p] + c) ^ 

(/ + T^E[sT, o\) n (e + T^E[sq,p] +b)^^ 

{d + TpE[r, o]) n (/ + T^E[sr, g] + c) ^ 



(e + T^V[sq]) n{d + T^F[sp, g, r] + a - 6) ^ 
(/ + T^.E[sr, o\) n (e + T^E[sq,p] + 6 - c) ^ 
[d + TpE[q, r]) n (/ + T^E[sr, o] + c - a) 7^ 
F[so, g, r] n (/ + T^l^[sr] + c - a) 7^ 

optimize (3 free parameters) 

The lattice vectors are 

a (+2.59649957, +0.85417034, +1.06467162) 
h (-0.14944338, +2.53196005, +0.58678219) 
c (-0.02900844, -0.61508232, +2.49547175) 

d (+0.02462012, -0.05626441, +0.08088453) 
e (-0.03576085, -0.12004584, -0.08428499) 
/ (-0.04295985, +0.03620318, -0.00675668) 

a + .07518716 
/? - .08591647 
7 + .04543238 

The lattice volume and packing density are 

y = det[a,6,c]17.82301085 
= 5- 74809657 

The clusters arc 

+F2 = +(5[o,p,g,r] UB[p,g,r, s]) 
-F2 = -(^[o,P,g,r] U5[p,g,r,s]) 
-Bi = -B[o,p,q,r\ 
+Bi = +S[o,p,g,r] 

The cluster vertices are 

o = (+1,+1,+1) 

P= (+1,-1,-1) 
q = (-1,+1,-1) 
r = (-1,-1. +1) 

* \ 3' 3' 3/ 

The neighbors of F2 are 

+ a — h F2 + 6 — c F2 + C — a 
F2-a + 6 F2-fe + c F2 - c + a 

F2 + ci 

e + f-F^-a e + f-F^-b e + f-F^- 
e — B]^ + a e — B^ + 6 e — B^ + c 
/ + Bi + a / + B1 + 6 / + Bi + c 
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The linear incidence conditions are 

V\p] n {F[o,q,r] + a - 6) 7^ 
V[q]n {F[o,r,p] + b-c) ^0 
V[r] n lF[s,p,q\+c-a) ^ 
E[o,p] n (E[s,r] +d) ^ 

g, r] n (e + / - g, r] - a) ^ 
F[s,r,p] n(e + /-F[s,r,p] -6) 7^0 
F[s,p,?]n(e + /-F[s,p,g]-c)^0 

F[o,p,g]n(e-F[o,p,g]+a) ^0 
F[o, q, r] n (e - F[o, g, r] + 6) ^ 
F[o, r, p] n (e - F[o, r,p]+c) ^0 

E[o,p] n (/ + E[g,r] +rf + a) ^ 
E[o,g]n(/ + E[r,p] + ci + 6) 7^0 
E[o,r] n if + E[p,q]+d + c) 

(/ + F[o, g, r] + d + a) n (e - F[o, g, r] + c) 7^ 

optimize (1 free parameter) 

The lattice vectors are 

a = ^ (-3069 + 11V396129, -2382 + 2^396129, -3369 - 1V396129) 
b = ^ (+17694 - 34^396129, +27651 - 37^396129, +3777 - 7^396129) 
c = ^ (-14625 + 23V396129, -25269 + 35^396129, -408 + 8^396129) 

= ^ (+22629 - 19V396129, -30903 + 65^396129, +4152 + 8^396129) 
e = ^ (+4131 - 5V396129, +2382 - 2V396129, -3777 + 7V396129> 
/ = (-2838 + 2V396129, +2865 - 7^396129, +1116 - 4V396129> 

The lattice volume and packing density are 

V = det[d + a,d + b,d + c] = (97802 181 - 132043^396129) 
= 6 • l/y = 11228544/(97802181 - 132043V396129) 

DIO The clusters are 

+E5 = +(BiUEj 

Bi B[t,u,v,w] 

E4 = B[o, p, V, w] U B[p, q, V, w] U B[q, r, v, w] U B[r, s, v, w] 
The cluster vertices arc 

o = V6 i-l-l-l) 
p = VQ{-l-l+l) 
q = VQ (0,0, +1) 

r^Ve i+l+l+l) 
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(+h+L-2) 
u = (-1.-1.-2) 
^; = (-l,+l,0) 
(+1,-1,0) 

The neighbors of Eg are 

E^ + a + b E5 + 6 + C Eg + c + a 
Eg-a-fe E5-6-C E5 - c- a 

d-Eg + a d-Eg + ft d-Eg + c 
d-Eg-a d-Eg-6 d-Eg-c 
d-E^ + a + b + c 
d — — a — b — c 

The linear incidence conditions are 

V[t]n{F[p,q,v] + a + b) ^0 
£^[r, v] n w] + 6 + c) 7^ 
r, w] n + c + a) 7^ 

F[t, u, v]n{d- F[t, M, + 6) ^ 

F[t, u, w]n{d- F[t, u,w]-c) 

F[o, p, v]n{d- F[o, p, v] - a) 7^ 

F[r, s, w]n{d- F[r, s, w] + a) 7^ 

F[q, r, v]r\{d- F[q, r, t>] + c) 7^ 

F\p, q, w]n{d- F\p, q,w]-b)^^ 

F[r, s, v]r){d- F[r, s,v] + a + b + c) $ 

F[o,p,w] n {d-F[o,p,w] - a-b-c)j^$ 

optimize (1 free parameter) 

The lattice vectors are 

a - ^ (+16713 + 14822\/6, -25455 + 13840v^, O) 
6 = ^ (-2656 + 561 V6, +14440 - 2385^6, -7285 - 4835v^> 
c = ^ (-2656 + 561v^, +14440 - 2385v^, +7285 + 4835v^> 
d = ^ (0, 0, -4339 + 1889v^> 

The lattice volume and packing density are 

y = 2 det[a, b, c] = 370146232 (''885808912 + 1639809563^6) 
= 10 • l/y = 29611698560/(23657426736 + 4919428689v^) 



